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Abstract. We report on new aspects of the star-forming region S235AB revealed through high-resolution observations at 
radio and mid-infrared wavelengths. Using the Very Large Array, we carried out sensitive observations of S235AB in the cm 
continuum (6, 3.6, 1.3, and 0.7) and in the 22 GHz water maser line. These were complemented with Spitzer Space Telescope 
Infrared Array Camera archive data to clarify the correspondence between radio and IR sources. We made also use of newly 
presented data from the Medicina water maser patrol, started in 1987, to study the variability of the water masers found in the 
region. S235A is a classical HII region whose structure is now well resolved. To the south, no radio continuum emission is 
detected either from the compact molecular core or from the jet-like structure observed at 3.3 mm, suggesting emission from 
dust in both cases. We find two new compact radio continuum sources (VLA-1 and VLA-2) and three separate maser spots. 
VLA-1 coincides with one of the maser spots and with a previously identified IR source (Ml). VLA-2 lies towards S235B and 
represents the first radio detection from this peculiar nebula that may represent an ionized wind from a more evolved star. The 
two other maser spots coincide with an elongated structure previously observed within the molecular core in the C^'^S line. This 
structure is perpendicular to a bipolar molecular outflow observed in HCO^(l-O) and may trace the associated equatorial disk. 
The Spitzer images reveal a red object towards the molecular core. This is the most viable candidate for the embedded source 
originating the outflow and maser phenomenology. The picture emerging from these and previous data shows the extreme 
complexity of a small (< 0.5 pc) star-forming region where widely diff'erent stages of stellar evolution are present. 
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1. Introduction 

This paper continues the study of the star-forming complex 
S235A-B (see Felli et al. 2004 and references therein), fo- 
cussing in particular on a deeply embedded Young Stellar 
Object (YSO) found between S235A and S235B, close to a 
water maser. The presence of the YSO is implied from the typ- 
ical signposts of early stellar evolution, including two molecu- 
lar outflows, a hot molecular core, a sub-millimeter peak, and 
a water maser. The YSO represents the youngest object in the 
star-forming complex. 

The site morphology derived from all previous observations 
is summarised in Fig. [I] S235A is a small optical nebulos- 
ity coinciding with a compact, but well-resolved, HII region. 
It appears as a less-evolved region, with respect to the more 
extended and diflTuse HII region S235 located further north, 
and probably is unrelated to the S235A-B complex. S235A lies 
at the northern edge of a molecular clump that represents the 
brightest peak of a more extended molecular cloud (Evans & 
Blair 1981; Nakano & Yoshida 1986; Cesaroni et al. 1999). 

Lying ~ W south of S235A (0.35 pc at the assumed dis- 
tance of 1.8 kpc; Nakano & Yoshida 1986), at the SW edge 
of the molecular core, S235B is a smaller diflTuse nebulosity 
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detected both in the optical and near-IR, exhibiting a near-IR 
excess and intense emission in optical and IR hydrogen lines 
(Ha, Bry), but without a radio continuum counterpart, making 
it a rather peculiar object. In Ha, it consists of an unresolved 
peak superimposed on a circular nebula that is ~ 10'' in diame- 
ter (Krassner et al. 1982; Alvarez et al. 2004). S235B appears to 
be a young star with an expanding ionized envelope surrounded 
by a diflTuse nebulosity (Felli et al. 1997). 

A large-scale molecular outflow was first found in ^^C0(1- 
0) by Nakano & Yoshida (1986) with a resolution of ~ 14''; it 
was centred at S235B and aligned in a NE-SW direction, about 
35" (0.3 pc) in length. Felli et al. (1997) confirmed the blue 
lobe of the outflow in ^^C0(2-l) with a resolution of ~ 11", 
but failed to detect the red lobe. 

Near-IR images revealed a highly obscured stellar cluster 
between S235A and S235B, i.e. centred on the water maser, 
with several sources with IR excess (Felli et al. 1997), in par- 
ticular source Ml, which exhibits the largest near-IR excess. 
This source has been suggested to be the candidate YSO sup- 
plying energy to the -60 km s"^ water maser, but with great 
uncertainty since it lies more than 5" to the south (see also 
Fig.ID. 

The above picture was clarified by Felli et al. (2004), 
who presented high-resolution (between 2" and 4") mm line 
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Fig. 1. Overview of the S235A-B star forming complex before the observations presented in this paper. Left: overlay of the 
single-dish HCN map (contours) from Cesaroni et al. (1999) with the ^-band image of the S235A-B region from Felli et al. 
(1997). Right: overlay of the interferometric maps of the 3.3 mm continuum emission (thick contours and grey scale) with the 
HCO+(1-0) outflows (thin solid line: V < -20.9 km dashed line: V > -15.5 km s"^) from Felli et al. (2004). In both, a 
cross marks the location of the H2O maser at -61.2 km s"^ from Tofani et al. (1995). The infrared source with the largest near-IR 
excess detected in previous works, Ml, is also indicated. The offsets between Ml, the water maser, the HCN peak, and the mm 
continuum peak are all larger than the position uncertainties. 



(HCO+, C^^s, H2CS, SO2, and CH3CN) and continuum obser- 
vations, together with far-IR observations. 

A compact molecular core (hereafter, the mm core) was 
found both in the mm continuum (hot dust emission, r^ust = 
T'cHgCN ~ 30 K) and in the molecular lines, peaking close 
to the water maser position and well- separated from S235A 
and S235B. Two molecular outflows were found in HC0^(1- 
0) centred on the mm core. One of them (hereafter, the NE- 
SW outflow) is aligned along the same NE-SW direction of 
the large-scale outflow detected by Nakano & Yoshida (1986). 
It spans ~ 0.4 pc, and has an estimated mass of 9 M© and 
a mechanical luminosity > 19 L©. The other (hereafter, the 
NNW-SSE outflow) is more compact and aligned in a NNW- 
SSE direction. It spans ~ 0.3 pc and has an estimated mass of 
4 M0 and a mechanical luminosity of ~ 0.3 L©. 

Felli et al. (2004) derived an upper limit of ~ 10^ L© for 
the bolometric luminosity of the mm core and suggested the 
presence of an embedded intermediate-mass YSO driving the 
NE-SW outflow and supplying the energy for the -60 km s"^ 
water maser. Studying the high velocity red and blue emission 
of C^^S(5-4) towards the mm core, they also found a compact 
structure with a velocity gradient perpendicular to the NE-SW 
outflow that might represent the signature of a circumstellar 
disk. 

An elongated structure (called a "jet") protruding from the 
mm core and coinciding with the blue lobe of the NNW-SSE 
outflow was also detected in the continuum at 3.3 mm (see 
Fig. [1] right) and, with a lower signal-to-noise ratio, also at 
1.2 mm. The spectral index of the jet, a ~ 0.6 (defined as 
Sy oc y^) is rather uncertain, but diff'erent from that of the mm 
core, a ~ 2.5, suggesting that the emission might arise in an 



ionized wind rather than being due to dust. Radio jets are often 
observed at the base of molecular outflows (Rodriguez 1997; 
Anglada et al. 1998; Beltran et al. 2001) in both low-mass and 
high-mass star-forming regions (e. g., Rodriguez 1996). They 
are characterized by spectral indices in the range -0.1 to ~ 1 
and are elongated in the outflow direction. Expanding ionized 
envelopes also have a spectral index a = 0.6 (Panagia & Felli 
1975). 

Water masers are one of the most reliable signposts of early 
phases in star formation (see e.g. Tofani et al. 1995) since 
they provide the best indication of the position of the required 
powering source, i.e. the YSO. They occur both in low-mass 
(see, e.g., Furuya et al. 2001, 2003) and in high-mass (see 
Church well 2002, and references therein) star-forming regions 
and are often found to be associated with outflowing matter. 
Sometimes they are also found in close association with radio 
jets (e.g. Gomez et al. 1995). 

The presence of a water maser in this region had been 
known since the observations of Henkel et al. (1986) and 
Comoretto et al. (1990), but with insufficient spatial resolu- 
tion to properly locate it in the region. Only with Very Large 
Array (VLA) cm line observations (Tofani et al. 1995) and in- 
terferometric mm observations of the continuum sources (Felli 
et al. 2004), was the location of the water maser in-between 
S235A and S235B, almost coincident with the mm core, firmly 
established. This proved that the water maser is not associated 
with either of the two nebulosities and that a local early type 
star, presumably the YSO within the mm core, is needed for 
its excitation. In the VLA observations, maser emission was 
only searched for in a limited velocity range, around -60 km 
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s"\ since at the time this was the only component detected by 
single-dish observations. 

The water maser in S235A-B has been monitored with the 
Medicina radio telescope since 1987, with coverage ~ 4 times 
per year since 1993. These observations revealed at least three 

separate velocity components: the one already known at 60 

km s"\ one between -20 to -30 km s"\ and one between -10 
to 10 km s"^. The last two are always very weak (max 10-20 
Jy) and exhibit strong variations. Whether they are all related to 
the same YSO or to separate ones could not be established from 
single-dish observations because of low spatial resolution. This 
drove us to carry out new VLA line observations at 22 GHz 
covering the whole maser velocity interval. All three velocity 
components found in the Medicina data were active at the time 
of this VLA observation. 

The S235A-B region contains other masers, namely 
methanol (CH3OH) and SiO (Nakano & Yoshida 1986; 
Haschick et al. 1990; Harju et al. 1998). Kurtz et al. (2004) 
included S235A-B in their recent VLA survey of the CH3OH 
maser line at 44 GHz. This is a class I methanol source, and 
it is believed to trace outflow activity. These authors found a 
cluster of 6 CH3OH masers spread over an area of a few square 
arcsec around the water maser spot at 60 km s"^ . 

To clarify the nature of the embedded YSO and its rela- 
tion with the outflow found in the S235A-B region, we have 
performed an extensive observational program using the VLA 
and the Medicina radio telescopes, complemented by archival 
Spitzer data. 

The two primary goals of the new VLA continuum observa- 
tions were: 1) to clarify the nature of the "jet" and to derive its 
spectrum over a larger frequency interval and 2) to search fur- 
ther for cm emission from ionized hydrogen in the mm core. 
VLA observations in the water maser line together with the 
Medicina patrol can give indications on the location and ac- 
tivity of the maser in the star-forming region. Finally, for a 
better understanding of the precise correspondence between IR 
and radio sources, in particular the precise role of Ml and the 
possibility of detecting IR emission from the mm core, archive 
Spitzer-IRAC observations of the S235A-B region in the four 
wavelengths (3.6, 4.5, 5.8, and 8 jim) were retrieved and ana- 
lyzed. 

In Sect. 121 we describe the observations, and in Sect. Owe 
present the results, while in Sect. we discuss our findings 
and how they enrich our current understanding of the S235A-B 
region. In Sect.|5l the main results are summarised. The reader 
can refer to Fig.[T6|for a comprehensive sketch of the S235A-B 
star forming-region, including the latest data. 

2. Observations and data reduction 

2.1. VLA observations 

The observations carried out with the VLA of the National 
Radio Astronomy Observatory (NRAO)^ were made with the 
VLA in the C configuration on February 26 and on March 7, 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



2004. The location of the water maser at -60 km s ^ was used 
as phase centre; its coordinates are: 
a(2000) = 05^0^^53.42^ ^(2000) = 35°4r48r8. 
The observations consisted of calibrator-target-calibrator 
scans; in the Q-band (0.7 cm) the fast- switching mode was 
used. The phase calibrator was 0555+398, with boot-strapped 
fluxes of 5.3 Jy (6 cm), 4.55 Jy (3.6 cm), 2.71 Jy (1.3 cm), and 
1.85 Jy (0.7 cm); the absolute amplitude calibrator was 3C147 
(0542+498), with assumed fluxes of 7.9 Jy (6 cm), 4.8 Jy (3.6 
cm), 1.8 Jy (1.3 cm), and 0.9 Jy (0.7 cm). The pointing accu- 
racy was checked on 0555+398 (at 3.6 cm) every hour. Line 
observations at 1.3 cm were carried out with the correlator set 
to 64 channels and a bandwidth of 6.25 MHz. This resulted in 
a velocity resolution of 1.3 km s"^ over a range of 84 km s"\ 
sufficient to cover the whole velocity range of the water masers. 
The receivers were tuned so as to have the band centred at -30 
km s"^ with respect to the Local Standard of Rest. A summary 
of the observations is contained in Table [I] 

Since the forthcoming analysis is based on a comparison 
of the position of the radio sources with those of near- and 
mid-IR sources, we checked whether the radio coordinate sys- 
tem is consistent with the near-IR coordinate system. To this 
end, we searched the 2MASS point source catalogue for a near- 
IR counterpart of our phase calibrator, finding an object at the 
same position as the calibrator within O'.' 1 . This gives us confi- 
dence that the two coordinate systems are consistent with each 
other within this limit. 

2.2. Medicina observations 

The single-dish Medicina radio telescope^ (HPBW If 9) obser- 
vations are part of a monitoring project of a large sample of 
water masers. An autocorrelator with 1024 channels and a 10 
MHz bandwidth is usually employed. The typical sensitivity 
for a 5 -minute integration is of the order of 1 Jy, and the av- 
erage calibration error is about 20%. For a more detailed de- 
scription of the radio telescope and the relevant parameters of 
the water masers patrolling, we refer to Valdettaro et al. (2002) 
and Brand et al. (2003). 

2.3. Spitzer-IRAC observations 

Spitzer-IRAC observations of a large area around the S235A-B 
region in the four wavelengths (3.6, 4.5, 5.8, and 8 jim) were 
extracted from the Spitzer public archive. The observations are 
part of the GTO program 201 "The Role of Photodissociation 
Regions in High Mass Star Formation" (Principal Investigator 
G. Fazio). Integration time is 12s at all filters. 

The positional accuracy is better than V\ Point-source 
FWHM resolutions range from ~ V!6 at 3.6 jim to ~ V!9 at 8.0 
yum. The IRAC bands are large and may contain various fea- 
tures, depending on the environment being observed. Among 
these, the most important for our case are polycyclic aromatic 
hydrocarbon (PAH) features (3.6, 5.8, and 8.0 bands) and the 
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Radioastronomy Institute of INAF, Italy 
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Table 1. Summary of VLA observations. 



Date 


Frequency 


Synthesized beam^ 


Largest 


Int. time 


rms 


Notes 






Size 


PA 


Angular Scale 










(GHz) 


(arcsec x arcsec) 


(degree) 


(arcsec) 


(sec) 


(mJy/beam) 




07/03/2004, 26/02/2004 


4.75 


4.8x4.3 


-8.2 


300 


1200 


0.09 




07/03/2004 


8.45 


2.8 X 2.6 


-20.8 


180 


1800 


0.08 




26/02/2004 


23 


1 X 1 


-40 


60 


7224 


0.03 


cont. 


07/03/2004 


23 


1 xO.9 


-83.7 


60 


434 


0.05 


line 


07/03/2004 


45 


0.6 X 0.6 


88.6 


43 


17590^ 


0.08 





^ natural weighting 
^ fast switching 



Bm line (4.5 band). An overview of IRAC is given by Fazio et 
al. (2004a). 



3. Results 

3.1. Spitzer-IRAC observations 

A colour coded image of the region covering S235A-B, ob- 
tained by combining 4.5, 5.8 and 8.0 jim observations, is shown 
in Fig.|3 

The two diffuse nebulosities S235A and S235B clearly 
dominate the extended emission. 

Most of the point sources detected in the ^-band and shown 
in Fig.Qare also present here. Two sources deserve more atten- 
tion in view of the present work: 1) Ml, which is detected at all 
bands, and 2) a source (hereafter, S235AB-MIR) detected only 
at 4.5, 5.8, and 8.0 yum, which is present in the area of the mm 
core. Both sources are marked in Figs.|2|and O where we show 
an overlay of the mm core (at 1.2 mm) with the 8.0 jim IRAC 
image. S235AB-MIR Hes Vf5 to the south of the peak of the 
mm core. 

We performed aperture photometry on the IRAC sources 
found near to the mm-core by using DAOPHOT in IRAK For 
all four bands, we selected a radius of 2 pix (~ 1 FWHM; 1 
pix ~ V!2) and a 2-pix wide annulus with an inner radius of 4 
pix, to account for the highly variable background. We applied 
aperture corrections as estimated from the IRAC PSFs retrieved 
from the Spitzer Web Page (http://ssc.spitzer.caltech.edu/obs/). 
To derive the photometric zero points, we used the zero- 
magnitude fluxes given by Fazio et al. (2004b). 

Photometry was done only on the most relevant sources 
present in the area of the S235A-B cluster, which are indicated 
in Fig. |2l Positions and flux densities in the four IRAC bands 
of the sources labelled in Fig.|2|are given in Table|2| 

Colour-colour plots are shown in Fig. |4l In the [5.8]- 
[8.0]/[3.6]-[4.5] plot, the region occupied by Class I and Class 
II objects (Allen et al. 2004) is indicated. Ml has colours typi- 
cal of a Class I object, while S235AB-MIR is located in a part 
of the plot redward of that occupied by YSOs of mass ~ 6 M©, 
according to Whitney et al. (2004), i.e. the region correspond- 
ing to [5.8]-[8.0] > 1 and [3.6]-[4.5] > 1. 




Fig. 3. The 1 .2 mm core (contours) from Felli et al. (2004) over- 
laid on the 5.8 jim Spitzer-IRAC image (grey scale, S235A and 
S235B have been saturated to show the weak emission from 
S235AB-MIR). The positions of the mm core, S235AB-MIR, 
and Ml are indicated. 

3.2. Radio continuum 
3.2.1. S235A 

The VLA radio continuum maps are dominated by the emis- 
sion from the compact HII region S235A. We calculated the 
integrated flux at all frequencies over the area within the 3cr 
contour of the emission at 6 cm (i.e. at the lowest resolution); 
the obtained values are listed in Table [3l The integrated flux at 
6 cm is in very good agreement with previous measurements, 
yielding 270 ± 40 mJy. Israel & FelH (1978) found 220 ± 30 
mJy at 21 cm and suggested a partially thick emission at 6 cm, 
but the ratio of fluxes at 6 and 3.6 cm that we derive from our 
data is in agreement with an optically thin emission. The pres- 
ence of an ionizing star of spectral type BO. 5 derived from the 
radio fluxes in previous works (e.g. Felli et al. 1997) is there- 
fore confirmed. In Fig. 13 we show the maps of S235A at 6, 3.6, 
and 1.3 cm. The 3.6 cm map is overlaid on the 5.8 fim Spitzer- 
IRAC image and shows a well-resolved spherical shape. The 
Spitzer-IRAC image reveals the ionizing star S235A* at the 
centre of the nebula, previously detected in the K band (Felli 
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Fig. 2. Three colour (4.5 blue, 5.8 green, and 8.0 yum red) image of the S235A-B region. The labels indicate the position of Ml, 
S235AB-MIR, and the other stars of the cluster for which we performed photometry. 



et al. 1997). The main feature of the radio maps is the asym- 
metry of the isophotal contours in the SE-NW direction. In the 
IRAC images, the morphology clearly indicates the presence of 
a brighter ridge SE of S235A*. 

The contours at 1.3 cm outline the brightest parts of the 
radio ridge. In the original maps, the emission is more frag- 
mented because of the high resolution and low surface bright- 
ness. In Fig. 13 the map has been smoothed to a resolution of 
3'\ At 1.3 and 0.7 cm, the radio fluxes are lower than expected 
from an optically thin emission, as was also found by Felli et 
al. (2004) at 3.3 mm. We attribute this to an instrumental efl'ect 
caused by the filtering of extended structures in the interfero- 
metric observations (see Table[l]). 

3.2.2. The mm core, the jet, and the radio compact 
sources VLA-1 and VLA-2 

At none of the 4 VLA wavelengths were we able to detect emis- 
sion from the mm core where the presence of an intermediate- 
mass YSO is suggested by the mm observations. This excludes 
any thermal emission from a UCHII region associated with the 



YSO with a flux density above the noise level given in Table (3] 
and implies that emission from the core is dominated by dust. 
At the same time, our non-detection is not in contradiction with 
the existence of a dust core. In fact, extrapolating at 0.7 cm, the 
flux density measured at 3.3 mm (20 mJy) with a spectral index 
of 2.5 (Felli et al. 2004), we obtain 3 mJy for the flux expected 
from the core. Felli et al. (2004) estimate that the core diameter 
in the continuum at the highest resolution (1.2 mm) is ~ 3''; 
hence, assuming that all the emission is uniformly distributed 
in a circle of V!5 in radius and using the synthesized beam size 
at 0.7 cm given in Table[l] we obtain an expected flux of ~ 0.12 
mJy/beam, i. e. < 2cr (see Table[T]). Hence, dust emission at 0.7 
cm could be present below the sensitivity limit of our observa- 
tions. 

Similarly, no extended emission from the jet was observed 
at any of our four frequencies. While this could be an efl'ect 
of over-resolution and low surface brightness of the jet at the 
shortest wavelengths, it definitely rules out the hypothesis of an 
ionized jet at 6 and 3.6 cm. In fact, at 6 cm, the flux density per 
beam area extrapolated from the 3.3 mm flux in the hypothesis 
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Table 2. MIR fluxes for the Spitzer-IRAC sources towards the mm core. 
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2.7 + 0.6 


< 20 


<59 


13 (S235AB-MIR) 


05M0'^53.4^ 


35041/49- 
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5.0 + 0.3 


14 + 1 
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14 
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35°4r53" 
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15 
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^ N-E of a close by source only visible in the 5.8 \m\ band. 



of an ionized jet (i.e. using a spectral index a = 0.6) would be 
a factor of 3 higher than the upper limit quoted in Table [3l 

At 1.3 and 0.7 cm, where the resolution is higher, two 
nearly unresolved sources are present. They have been named 
VLA-1 and VLA-2 and are indicated in Fig.|6| where we show 
the VLA map at 1.3 cm (contours) overlaid on the Plateau de 
Bure map at 1.2 mm (grey scale) from Felli et al. (2004). We 
derived the integrated fluxes of the two sources within the 3(T 
level on the maps obtained with natural weighting. At 3.6 and 6 
cm, the two sources fall within the sidelobes of S235A. This re- 
sults in a noise higher than that predicted from the total integra- 
tion time (see Table [I]) . Using difl'erent weightings to partially 
filter out the extended emission does not yield any significant 
improvement in the measurements. The flux densities are listed 
in TableO the upper limits at 6 cm (for both) and at 3.6 cm (for 
VLA-2) refer to a point source. 

VLA-1 lies ~ W south of the mm core and is located 
within the jet, almost along its axis (see Figs.[l]and|6|). It co- 
incides with a small component of the fragmented jet observed 
at 1.2 mm and, most noticeably, with the K-band source Ml. 

VLA-2 lies within the boundary of the S235B nebulosity. 
VLA-2 represents the first radio detection of this peculiar re- 
gion. Our flux densities of ~ 0.5 mJy (see Table [IJ are not in 
conflict with the previously derived upper limits of 5 mJy at 6 
cm (10'' beam; Israel & FelH 1978) and 0.3 mJy at 3.6 cm (a.'l 
beam; Tofani et al. 1995). 

We have checked the probability that VLA-1 and VLA-2 
are background sources. The expected number of extragalactic 
sources at 1.3 cm in the field of view of Fig. |6|(~ 1 square 
arcmin) based on Eq. (All) of Anglada et al. (1998) is ~ 
0.05, making this possibility very unlikely. 

3.3. H2O masers 
3.3.1. VLA observations 

In the selected velocity range (roughly from -70 to 10 km 
s"^), three maser spots were detected above the 5cr noise. We 



have determined the flux densities and positions of the maser 
spots by 2-dimensional Gaussian fits in each channel. All ve- 
locity components in the same spot are spatially unresolved. 
Coordinates of the three maser spots and flux density for each 
velocity peak are listed in Table |4| their relations to the other 
features present in the area are shown in Figs. l6l and fT3l 

It is important to note that the three maser spots cover dif- 
ferent velocity ranges, as shown in Fig. Q so that there is no 
velocity overlap among the three spatial components. 

One of them (S235A-B-H20/3) coincides, within ~ a.'5, 
with that found with the VLA at ~ -60 km s"^ by Tofani 
et al. (1995) and emits in the same velocity range. The other 
two (S235A-B-H20/1 and S235A-B-H20/2) occur at radial 
velocities that had not been searched for in the previous VLA 
observation because at that time they were not detectable in 
single-dish observations, but they have since been revealed in 
the Medicina patrol. 

The maser luminosity for each spot was obtained by inte- 
grating the line emission within the respective velocity ranges. 
The results are listed in Table |5j along with the corresponding 
velocity range. The H2O luminosity is typical of masers as- 
sociated with far-Infrared (FIR) sources of 10^ - 10^ Lq (see 
Palagi et al. 1993). This is consistent with the upper limit of the 
bolometric luminosity inferred for the mm core by Felli et al. 
(2004). 

From Figs. |6| and El S235A-B-H20/2 clearly coincides 
with VLA-1 and Ml and does not seem to be directly related 
to the mm core. The other two spots (S235A-B-H20/1 and 
S235A-B-H20/3) are very close to the mm core and are per- 
pendicularly aligned to the NE-SW outflow, ~ apart from 
each other. The possibility that they might be tracing a disk or 
torus perpendicular to the NE-SW outflow will be examined in 
Sect. El 
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Table 3. Radio continuum fluxes. Upper limits are 3cr noise 
levels. 



Source 


Position 




^y(mJy) 






Qf(2000) 


(5(2000) 


0.7 cm 


1.3 cm 3.6 cm 


6 cm 


S235A 


05M0'"52.70^ 


35°42'21" 


114 + 3^ 


172 + 1 248 + 2 


257 + 7 


VLA-1 


05M0"^53.60^ 


35°41'43" 


0.59 + 0.08 


0.44 + 0.03 0.39 + 0.08 


<0.27 


VLA-2 


05M0'"52.40^ 


35°4r30" 


0.48 + 0.08 


0.47 + 0.03 < 0.24 


<0.27 



Corrected for primary beam attenuation. 



D 




10' 



^50 



Vlsr (^^^^/s) 



10' 



50 



Fig. 9. Time-velocity-intensity plot of the water maser emission from S235A-B observed with the Medicina radio telescope. The 
starting date is March 31, 1987. The dates of the first VLA observation by Tofani et al. (1995) and of the present observations are 
indicated (month/year) with an arrow and their velocity ranges are enclosed within a long rectangle. The three maser spots are 
indicated by bracketing the velocity ranges with vertical dashed lines. The vertical solid line defines the velocity of the molecular 
cloud. The black areas are time- velocity regions with no observations. 



3.3.2. Medicina observations 

The diff'erent velocity ranges of the maser spots detected with 
the VLA make the low-resolution Medicina observations a use- 
ful tool to follow the evolution of each spot, with no confu- 
sion arising from the fact that they are all within the Medicina 
beam. The upper envelope of the water maser emission using 
all the available single-dish observations is shown in Fig. [Sl 
When compared with Fig.Cl it shows that the emissions from 
the three spatially separated VLA maser spots occur at very dif- 
ferent velocity ranges, being thus separated (although spatially 
unresolved) in the Medicina observations, as well. 

In Fig.|9| we show the time- velocity-intensity plot from the 
Medicina patrol. An indicative value of the noise level in these 
observations throughout the entire period is of the order of 1-2 



Jy. The starting date corresponds to March 31, 1987. The patrol 
is sparser at the beginning. After 1992, there are about 4 obser- 
vations every year. Following the separation in velocity of the 

three maser spots, each velocity component (at 60, 25, 

and ~ km s"^) is labelled in Fig. |9| with the corresponding 
VLA name. The dates of the two VLA observations are indi- 
cated. The systemic velocity of the molecular cloud (-17 km 
s"^) is also indicated. Only S235A-B-H20/2 occurs at a ve- 
locity close to that of the thermal molecular lines; the other 
two components (S235A-B-H20/1 and S235A-B-H20/3) emit 
well outside the width of the thermal molecular lines. 

The biggest flare occurred from the component S235A-B- 
H20/3 at the time of the first VLA observation. The source then 
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Vlsr (knn/s) 

Fig. 10. Time- velocity-intensity plot of the water maser emission from S235A-B observed with the Medicina radio telescope 
over the velocity range from -10 to +10 km s"\ which corresponds to the maser spot S235A-B-H20/1. 



Table 4. Water masers: fluxes and positions. 



Name Position 



S235A-B 


Qf(2000) 


(5(2000) 


^LSR 

(km s-^ 


(Jy) 


H20/1 


05M0'"53.27^ 


35°41'5O/0 


7 


1.20 








3 


0.50 


H20/2 


05M0'"53.63^ 


35°4r43r0 


-18 


0.16 








-25 


0.52 








-29 


1.66 


H20/3 


05M0"^53.38^ 


35°4r48r6 


-58 


0.37 








-62 


0.23 








-64 


0.22 








-68 


1.44 



Table 5. Water masers: line luminosity. 



Name 


T ci 
^H20 


AVlsr 




(10-^ Lo) 


(km s 




S235A-B-H20/1 


3.1 


9to- 


-1 


S235A-B-H20/2 


3.3 


-17 to 


-30 


S235A-B-H20/3 


4.1 


-55 to 


-71 



^ All masers are assumed to be at the distance of the S235A-B com- 
plex. 

disappeared below the noise and came up again just shortly 
before the second VLA observation. 



Component S235A-B-H20/2 was undetectable for most of 
the time and appeared above the noise only after 2000. 

Component S235A-B-H20/1, although always rather 
weak, is present most of the time. Its most noteworthy aspect is 
that the velocity changes up to ± 5 km s"^ around a mean value 
of ~ km s"^ To better illustrate this effect, the time- velocity- 
intensity plot for the velocity range from -10 to 10 km s"^ is 
shown in Fig.fTOl 

To make sure that the change in velocity is a real eff'ect and 
not an instrumental one and to provide an independent estimate 
of the accuracy on the velocity, we have examined another wa- 
ter maser (G32.74-0.08), that is also included in the Medicina 
patrol, and that was chosen because it is characterized by a sin- 
gle, narrow, and intense velocity component. For this maser, 
the velocity of the peak displays a maximum deviation from 
the mean value over the entire period of < 0.1 km s"^ a factor 
of 50 smaller than the velocity spread observed in S235A-B- 
H20/1. 

4. Discussion 

4.1. S235A 

In Fig. 13 the overlay of the 3.6 cm map with the 5.8 jim Spitzer 
image shows that the peak of the IR emission occurs in a shell 
outside the boundary of the radio emission. This proves that 
PAH and thermal dust emissions are mostly located beyond the 
ionization front in the Photodissociation Region. 
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In Fig. [H] we show the overlay of the HCO^ integrated 
emission around -17 km s"^ with the 5.8 yum Spitzer-IRAC im- 
age. The overlay indicates that the lower contours of the molec- 
ular emission closely follow the outer boundary of S235 A, sug- 
gesting that the HII region and the molecular cloud are inter- 
acting. It is also worth noting that the mm core is just outside 
the S235A boundary. This situation resembles that of some 
well-known cometary- shaped UCHII regions such as G29.96- 

0. 02 and G34.26+0.15 (Reid & Ho 1985; Wood & Churchwell 
1989; Fey et al. 1995), which face a density peak of the molec- 
ular clump enshrouding them (Maxia et al. 2001; Gibb et al. 
2004; Watt & Mundy 1999). In a number of cases, it has been 
found that such a peak coincides with a hot molecular core, 
where massive star formation is going on (Cesaroni et al. 1998; 
Heaton et al. 1989; Garay & Rodriguez 1990). In our case, 
the temperature of the molecular core is ~ 30 K (see Felli et 
al. 2004), well below the typical temperature of hot molecular 
cores (see Kurtz et al. 2000), and the HII region S235A has 
an asymmetric rather than a cometary shape; nevertheless, the 
interaction between the ionized gas and the molecular cloud, 
as traced by the overall structure of the region, suggests that 
in this case also, one is observing an active burst of star forma- 
tion where different evolutionary phases (from cores to evolved 
HII regions) co-exist. Whether the star formation episode in 
the molecular core has been triggered by the expansion of the 
nearby HII region S235A remains an open issue. 

4.2. yLA-1 and the jet 

Our radio continuum observations were unable to detect the 
elongated structure that had previously been found at 3.3 mm, 
which also coincides with the blue lobe of the NNW-SSE 
HCO^(l-O) outflow. Instead, a compact radio source, VLA- 

1, was found coincident with a small blob of emission in the 
central part of the 1.2 mm "jet", close east) to the weaker 
molecular peak at -19 km s"\ called C19 (Felli et al. 2004). 
More importantly, VLA-1 coincides with the newly found wa- 
ter maser S235A-B-H20/1 and with Ml, as shown in Fig.fHl 
where the 1.3 cm VLA map is overlaid on the Spitzer-IRAC 3.6 
[im image. Ml, which had been previously assumed to be as- 
sociated with the -60 km s"^ water maser (Felli et al. 1997), is 
instead associated with a a radio continuum source and a new, 
separate water maser at 30 km s"^ . 

Although the errors on the flux densities are large, the val- 
ues given in Tabled for VLA-1 are consistent with a spectral 
index a ^ 0.2. This value is typical of partially optically thick 
free-free emission. 

Since the present morphology no longer supports an ion- 
ized jet interpretation, we have to consider the alternative pos- 
sibility that VLA-1 is an independent UCHII region, a hypoth- 
esis that should also be considered, in light of the precise as- 
sociation of the water maser with Ml. A lower limit to the to- 
tal number A^Ly of ionizing photons can be obtained from the 
1.3 cm radio flux (the one with the highest signal-to-noise ra- 
tio) by assuming that the HII regions are optically thin (Mezger 
1978). We obtain A^Ly ~ 10^^ s'^ , typical of B2-B3 ZAMS stars 
(Panagia 1973). 



What may have caused the apparent disagreement between 
the cm radio observations which postulate the presence of a 
UCHII region and the mm observation that had suggested a 
jet? Extrapolating the radio flux of VLA-1 according to 
(of ^ 0.2), we derive ~ 0.8 mJy at 1.2 mm and ~ 0.6 mJy at 3.3 
mm. The 1.2 mm flux is below Icr (~ 1 mJy) in the Plateau de 
Bure map, so it could not have been detected as a point source. 
The 1.2 mm flux of the "jet", 13 mJy, was integrated over a 
much larger area defined by the 3.3 mm map. At 3.3 mm, the 
extrapolated flux of VLA-1 is at a 2(T level, again barely de- 
tectable as a point source, and in any case difficult to see since 
it would lie within the elongated 3.3 mm structure (about 10'' 
in the elongated direction and unresolved in the perpendicular 
direction), with 7 mJy of integrated flux. 

These contradictory aspects can be reconciled if the emis- 
sion from the elongated 3.3 mm structure comes from dust, 
perhaps that associated with the NNW-SSE blue outflow lobe 
or with the C19 molecular peak. The spectral index a = 0.6 
between 3.3 and 1.2 mm could be the result of the simulta- 
neous presence of a dusty jet with a steep spectral index and 
a weak UCHII region showing up only at lower frequencies, 
where dust emission is negligible. 

Is VLA-l/Ml the driving source of the NNW-SSE outflow? 
Due to the overlap with the NE-SW outflow, the centre of the 
NNW-SSE outflow is ill-defined and any argument based on the 
position of VLA-1 with respect to the centre of the NNW-SSE 
outflow is inconclusive. However, we can exclude any associ- 
ation of VLA-1 with the NE-SW outflow, given its ofl'set with 
respect to its axis. 

Finally, we note that one of the methanol masers 
(CH3OH/4) Hes close (2" southwest) to VLA-1 and the wa- 
ter maser S235A-B-H20/2 (see Fig.fT3t and has velocities in a 
similar range (from -16 to -21 km s"^ see Kurtz et al. 2004), 
suggesting a common origin. The velocities of the water and 
methanol masers are very close to that of the molecular core, 
indicating that in this case, the component of the motions along 
the line of sight of the masers with respect to the molecular core 
is negligible. 

4.3. yLA-2 and S235B 

The overlay of Fig. fT2l shows that VLA-2 lies at the centre of 
the S235B bright diffuse nebula (saturated in all IRAC bands). 
The Spitzer-IRAC images, in particular those at longer wave- 
lengths, indicate that this nebula is composed of two very close 
components (see Fig.|2|). VLA-2 lies at the centre of the south- 
ern and more extended one. Overall, the size of the mid-IR 
nebula is about 10'', similar to that observed in Ha, so that the 
optical-IR morphology is more reminiscent of S235 A (resolved 
HII region) than that of VLA-1 (IR and radio unresolved). 

Besides being detected in Ha, strong Bry emission from an 
unresolved source (i.e. < 3'') coincident with the ^-band point 
source in S235B had been reported by Krassner et al. (1982) 
and Felli et al. (1997), with an integrated line flux of ^(Bry) 
= (2.0 ± 0.4) X 10"^^ erg s"^ cm"^. The expected radio flux 
density from an optically thin HII region at 6 cm would have 
been greater than 200 mJy, where the lower limit accounts for 



10 



M. Felli et al.: The S235A-B star-forming region 



the fact that the Bry flux was not corrected for extinction. This 
is at odds with lower limits to the radio flux found in all previ- 
ous radio continuum observations, as well as with the present 
detection. In the past, this forced abandon of the hypothesis of 
a classical HII region (unless extremely optically thick) sug- 
gested that the Bry emission originates from an ionized ex- 
panding envelope around an early type star (Felli et al. 1997), 
in which case the ratio of radio-to-IR line emission would be 
reduced by about two orders of magnitude (Simon et al. 1983). 

Our measured fluxes represent the first detection of the ra- 
dio continuum from an underlying unresolved source. The ra- 
dio flux density at 3.6 cm expected from a fully ionized en- 
velope using F(Bm)/F(Bry) -0.8 and Eq. (20) of Simon et 
al. (1983) is 0.46 mJy. While this agreement with the observed 
value might be fortuitous in view of the many unknown pa- 
rameters involved (velocity of the wind, correction for extinc- 
tion, etc.) it clearly points out that an ionized envelope around a 
lower luminosity,more-evolved star remains the best interpreta- 
tion. However, it must be noted that the observed spectral index 
is smaller than the expected 0.6 value. 

Our flux densities are fully consistent with the previous up- 
per limits at radio wavelengths and with the unresolved nature 
of the S235B. The implied mass loss using Eq. (14) of Felli & 
Panagia (198 1) is 4x 10"^ Mq yr"\ quite large and indicative of 
a luminous star. Future radio recombination line observations 
with sufficiently high sensitivity may permit measurement of 
the line width to determine the velocity of the ionized wind. 
There are no masers, molecular peaks, outflows, or mm peaks 
associated with S235B, all of which suggest that S235B-VLA- 
2 is more evolved than the mm core and VLA-1. 

The difl'use Ha and IR emission from S235B may be at- 
tributed to reflected light from ionized stellar envelopes. 

Finally, what is the luminosity and mass of the star embed- 
ded in S235B? We cannot use the Spitzer-IRAC observations 
because they are saturated. Instead, we used the J, H, and K 
magnitudes from Felli et al. (1997) and the MSX fluxes (see 
Egan et al. 1999). The integral of these values gives 410 Lq, 
which must be considered to be a lower limit because the FIR 
part of the spectrum is not taken into account in our calculation. 

4.4. The mm core and the S235A-B-H20/1 and 
S235A-B-H20/3 masers 

As was already pointed out, an important result of the cm ra- 
dio continuum observations is the lack of emission from the 
mm core. In Felli et al. (2004), the luminosity of the embed- 
ded YSO had been estimated to be 10^ Lq. If this comes from 
a ZAMS star of spectral type B3 or earlier, to make the ra- 
dio free-free emission undetectable, the radius of the associated 
HII region must be less than -6.5 AU due to confinement by a 
density larger than 7 x 10^ cm"^. In any case, the detection of 
pure thermal dust emission from the core indicates a very early 
evolutionary phase of the embedded YSO. 

The new result provided by the Spitzer-IRAC observations 
is the detection of S235AB-MIR very close to the mm peak. 
The positions of the mm core and S235AB-MIR diff'er by 
~ Vf5, which is slightly greater than the error on the relative 



positions, so that at this stage it cannot be firmly established if 
S235AB-MIR represents the YSO embedded in the mm core. 
Its position in the colour-colour plots of Fig. |4l definitely puts 
this source in the Class I category. 

It is possible to check whether the S235AB-MIR fluxes be- 
tween 3.6-8 Jim are consistent with the spectral energy distri- 
bution (SED) of a heavily embedded (proto-)star of 10^ Lq. 
Felli et al. (2004) show that the non-detection of such a source 
in the K band towards the mm core implies Ay ^ 37 mag, in 
agreement with the derived H2 column density. Using the ex- 
tinction law found by Indebetouw et al. (2005), we derived up- 
per limits for the intrinsic fluxes (or upper limits) of S235AB- 
MIR in the four IRAC bands. Of course, S235 AB-MIR being a 
Class I source, most of the emission at these wavelengths arises 
from circumstellar matter rather than from the (proto-)star pho- 
tosphere. For this reason, we compared our results to the SEDs 
models for intermediate-mass Class I sources (star plus disk 
and envelope, Whitney et al. 2003; Whitney et al. 2004). We 
found that, within the limits of the many parameters of their 
models, the MIR SED of S235AB-MIR is consistent with a 
central star later than B3. 

The two northernmost water masers, S235A-B-H20/1 and 
S235A-B-H20/3, are located close to the mm core and emit 
over velocity ranges quite difl'erent from those of the other wa- 
ter and methanol masers, as well as those of the molecular 

cloud ( 17 km s"^). Figure[T4| (adapted from Fig. 21 of Felli 

et al. 2004) shows the NE-SW HCO+ outflow with the con- 
tours of a perpendicular bipolar structure traced by the wings 
of C^4s(5_4). This was interpreted by FelH et al. (2004) as the 
signature of a rotating disk around the YSO driving the out- 
flow. Strikingly enough, the two water masers are aligned in the 
same direction and very close to the C^'^S(5-4) structures: the 
blue-shifted maser (S235A-B-H20/3 at ~ -60 km s"^) lying 
towards the blue C^'^S(5-4) lobe (integrated from -21 to -19 
km s"^) and the red-shifted maser (S235A-B-H20/1 at ~ 7 km 
s"^) lying towards the red C^'^S(5-4) lobe (from -16 to -14 km 
s"^). However, in both cases the velocity of the water masers 
are more red or blue shifted than the corresponding C^'^S(5-4) 
lobes. 

A simple calculation based on the maser velocities shows 
that they cannot simply belong to a disk in Keplerian rotation 
around the YSO. From the difl'erence between the most ex- 
treme maser velocities (7 and -68 km s"^), we can infer a lower 
limit to the rotation velocity of 37 km s"^ If one assumes the 
half-distance between the two maser spots (~ V or 1800 AU) 
to be the orbital radius, then the mass needed to maintain such a 
rotating disk should be > 2500 Mq, much larger than the mass 
of the molecular core. 

This fact proves that the water masers cannot trace rotation 
about an embedded YSO. One possibility is that they are in- 
stead the signature of the interaction with the C^^S disk of high 
velocity outflowing material from the YSO. The higher (red 
and blue) velocities of the water maser with respect to C^^S 
and their closer position to the mm core may suggest that the 
acceleration of the outflowing material occurs in the immediate 
surroundings of the YSO. 

Another possibility is that the C^^S emission is not trac- 
ing a disk, but rather an outflow, whose high-velocity com- 
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ponent would be seen in the H2O maser lines. This scenario 
is more consistent with the common belief that water masers 
are strictly associated with jets powering molecular outflows 
(Felli et al. 1992), as observed, e.g., in the massive protostar 
IRAS 20126+4104 (Moscadelli et al. 2000, 2005). If this hy- 
pothesis is correct, it remains to be established whether the 
bipolar outflow seen in the C^'^S and H2O lines would be the 
same as the NNW-SSE outflow or a distinct one oriented ap- 
proximately in the same direction, but originating from a dif- 
ferent YSO. 

Of the two hypotheses presented above, we believe that the 
outflow origin for the C^^S and H2O emission is the most likely, 
given the tight association between outflows and water masers. 
However, at present it is impossible to rule out the possibil- 
ity that the C^^S emission is coming from a disk, as in the 
case of the Keplerian disk in IRAS 20126+4104 (Cesaroni et 
al. 2005). Class I methanol masers are believed to be excited in 
jets, so that a priori the 7 mm CH3OH masers imaged by Kurtz 
et al. (2004) in S235 could be used to establish the direction of 
the outflow and hence choose between the two hypotheses. As 
shown in Fig.[21 five of the maser spots cluster around the mm 
core suggesting that for them, too, the main source of energy is 
the YSO within the mm core. In particular, CH3OH/2 is very 
close to S235A-B-H20/3. However their velocities lie within 
a narrow range (from -15.9 to -21.0 km s"^) and the spots do 
not show a clear bipolarity, either in velocity or in their distri- 
bution. It is therefore difficult to associate the methanol maser 
emission to any precise geometry. The only conclusion is that 
they are unlikely to have the same dynamical origin as the two 
water masers. 

4.5. Water maser variability 

Some indication of what occurs in the immediate surroundings 
of the YSO embedded within the mm core may come from the 
variability of the two associated water masers. 

The emission at ~ -60 km s"^ (S235A-B-H20/3) reached 
its maximum (up to ~ 1 10 Jy) in 1992-1993, lasting at most 2 
years, then disappearing for most of the Medicina patrol, and 
finally reappearing just before the VLA observations. 

Most noteworthy is the emission from S235A-B-H20/1 
(between -10 and 10 km s"^), not only because of its long 
lasting presence and high variability, but also because of its 
velocity shifts, +5 km s"^. Figure [21 shows the LSR velocities 
of S235 A-B-H20/1 (obtained by Gaussian fitting) in the period 
1989-2005. 

While it is possible that the velocity variations are simply 
due to the random flaring of spots at difl'erent velocities, we 
shall investigate more physically appealing interpretations. 

To check whether the efl'ect is due to a rotational modu- 
lation, we tried to fit all the observed velocities (the source is 
above the noise for about 60% of the time) with a sine function, 
as would be expected for a maser spot on a rotating disk viewed 
edge-on. No significant evidence was found. 

The velocities in Fig. [21 are all redshifted with respect to 
the systemic velocity, but they are not all at random and seem 
to fan out at late periods. We identify two groups of points (la- 



belled as 1 and 2 in Fig. fTSt . which exhibit linear drifts of ve- 
locity away from the systemic value. A linear fit (dashed lines 
in Fig.fTSt provides velocity drifts of 0.93 and 0.98 km s"^ yr"\ 
respectively. Velocity drifts of this amount have been observed 
in other water masers that are stronger and less variable in in- 
tensity (Brand et al. 2003). It is tempting to explain the veloc- 
ity drifts with shocked material that is accelerated from a mean 
velocity ~ km s"^ by mass outflow from a central YSO. The 
lifetime of the accelerated spots is ~ 1000 - 2000 days (-3-6 
yrs) and could be related to the duration of ejection events from 
the YSO. 

The maser outflows from S235A-B-H20/1 clearly deserve 
further, proper motion studies with VLBI techniques. 

5. Summary and conclusions 

We have presented new, more sensitive high-resolution VLA 
cm radio observations of the S235A-B region, as well as the 
results of the Medicina water maser patrol (started in 1987), 
and archive Spitzer-IRAC observations. Several new aspects of 
this star-forming region emerge; they are illustrated in Fig.fT6l 
and summarised in the following: 

1. The radio-IR morphology of the S235A HII region con- 
firms that it is a classical HII region, optically thin in the 
cm range. It appears to interact with the molecular cloud 
and may have induced the formation of a second genera- 
tion YSO in the mm core. 

2. No cm continuum emission is detected from the molecu- 
lar core discovered by Felli et al. (2004). The lack of ion- 
ized hydrogen emission suggests a very early evolution- 
ary phase for the intermediate-luminosity embedded YSO, 
much before the appearance of a UCHII region. We have 
found a new source, S235AB-MIR, detectable only at 5.8 
and 8.0 yum and close to the mm core, in the archival Spitzer 
images. Given its position in the colour-colour plot, this 
could be the mid IR counterpart of the embedded YSO. 

3. We have observed no extended cm continuum emission 
from the elongated jet-like structure detected at 3.3 mm, 
suggesting that the putative 3.3 mm "jet" is due to dust and 
not to an ionized jet. 

4. We found two compact radio- sources: VLA-1 and VLA-2. 
Their spectral index is suggestive of partially thick free-free 
emission. 

5. VLA-1 is located at the centre of the elongated 3.3 mm 
structure and coincides with the near-IR source Ml. It is 
close to the secondary molecular peak at -19 km s"^. We 
estimate that VLA-1 could be a UCHII region associated 
with a B2-B3 star. 

6. We have discovered a water maser (S235A-B-H20/2) at the 
same location of VLA-1. A methanol maser (CH3OH/4) 
lies close by at a similar velocity. 

7. VLA-2 is at the centre of S235B and represents the first 
radio continuum detection from this source. Comparison 
with the near-IR hydrogen lines confirms that both emis- 
sions came from an ionized envelope. 

8. Two water masers (S235A-B-H20/1 and S235A-B- 
H20/3), with very different velocities (~ 7 and 60 km 
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respectively), are located close to the mm core and are 
aligned parallel to a structure found in C^'^S and perpendic- 
ular to the NE-SW outflow. 
9. Our single-dish observations with the Medicina radio tele- 
scope do not spatially resolve the three water masers de- 
tected with the VLA, although do provide variability infor- 
mation since they do not overlap in velocity. 

10. A high degree of variability in the water maser emission 
was found in all cases. We found changes of the LSR ve- 
locity with respect to the systemic velocity, up to ~ 5 km 
s~^ for S235A-B-H20/1. A possible interpretation would 
be velocity drifts due to shocked gas accelerated by the flar- 
ing activity of the YSO. 

1 1 . The duration of the acceleration, of an order of 3 - 6 yrs, is 
similar to the lifetime of the emission. 

This paper, together with the preceding ones resulting from 
our long-term study, reveals the simultaneous presence in the 
S235A-B complex of widely diflferent evolutionary phases, 
from the well-developed HII region S235A, to the peculiar 
object S235B-VLA-2, to the UCHII region VLA-1, to the 
mm core which harbors a YSO in a very early stage. The 
comparison of high resolution multiwavelength observations, 
combined with long-term monitoring of time variable phe- 
nomena, provides unique information on the nature of young 
(proto) stars, shedding new light on their interaction with the 
placental environment. 
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Fig. 4. Colour-colour plots of the most relevant sources in the 
area of the S235A-B cluster. The identifying numbers and sym- 
bols are the same as those used in Table|2land Fig.|2l In the bot- 
tom box, the regions occupied by Class I and Class II sources 
are enclosed with dotted and dashed lines, respectively. When 
only an upper limit to the flux density could be estimated in 
one of the bands, the corresponding point in the plot is marked 
by an arrow. 
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Fig. 5. Maps of the continuum radio emission from S235A 
at 6, 3.6, and 1.3 cm (contours). Levels are from 2 to 14 
mJy/beam by 2 mJy/beam for 6 cm, from 1 to 5 mJy/beam 
by 0.5 mJy/beam for 3.6 cm, and 0.1 to 0.5 mJy/beam by 0.1 
mJy/beam for 1.3 cm. The 1.3 cm map has been smoothed to 
a lower resolution to increase the S/N ratio. The 3.6 cm map is 
overlaid with the 5.8 jum Spitzer-IRAC image (grey scale). 
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Fig. 6. Map of the continuum radio emission at 1.3 cm of the 

area around the water masers (contours) overlaid with the map 
of the continuum 1.2 mm emission (grey scale) from Felli et 
al. (2004). Levels are from 0.12 (~ 4o-) to 0.48 mJy/beam in 
steps of 0.12 mJy/beam. The locations of the three H2O maser 
spots detected in this work are indicated by "+". Ml and the 
centre of S235B (asterisk) are also indicated. The dashed lines 
define the directions of the two outflows. The dotted circle is 
the primary beam at 1.2 mm. The synthesized beam at 1.3 cm 
is drawn in the bottom left-hand corner. 
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Fig. 7. H2O VLA spectra averaged on a circle \" in radius and 
centred on the locations of each of the 3 maser spots. The name 
of the spot is indicated above the corresponding spectrum. The 
intensity scale is off'set by 2 and 4 Jy for clearness of the dis- 
play. The vertical line defines the velocity of the molecular 
cloud. 
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Fig. 8. Upper envelope of all water maser spectra observed 
with the Medicina radio telescope towards S235A-B, until July 
2005. The three corresponding maser spots detected with the 
VLA observations are indicated. The vertical line defines the 
velocity of the molecular cloud. 
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Fig. 11. The HCO^ cloud at -17 km s"^ (contours) overlaid 
on the 5.8 jim Spitzer-IRAC image (grey scale). Note that the 
lower contours of the molecular emission closely follow the 
boundary of the IR emission, suggesting that the S235A HII 
region is interacting with the molecular cloud. 
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Fig. 12. Overlay of the 1.3 cm VLA map (full contours) with 
the Spitzer-IRAC 3.6 jim image (grey scale) and the 1.2 mm 
continuum (dashed contours). 
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Fig. 13. Position of the three water masers (H20/1-3 +) and of 
the six methanol masers (CH30H/1-6 x) detected by Kurtz et 
al. (2004). The 1.2 mm core (triangle), VLA-1 (square), and 
Ml are also indicated. Note that the symbol used to mark the 
location of the 1.2 mm core is smaller than the errorbars. 




Fig. 14. HCO+(1-0) emission from the NE-SW outflow (thin 
contours) overlaid with the blue- and red-shifted emission of 
C^4s(5-4) (thick contours) from Felli et al. (2004). The 
symbols mark the location of the water maser spots H20/1 

(~ 7 km s-i) and H20/3 ( 60 km s'^). The background 

is the 8.0 /im Spitzer-IRAC image (grey scale) and shows the 
position of S235AB-MIR, coincident with the southern water 
maser (H20/3). 
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Fig. 15. Velocities of the peaks of S235A-B-H20/1 from 
Gaussian fits to the data from the Medicina radio telescope. We 
have tentatively outlined, with dashed lines, two components 
(labelled as 1 and 2) whose velocities might be drifting during 
the period of our observations. 
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Fig. 16. Sketch (not to scale) of the star-forming region S235A- 
B in light of the new data presented in this paper. New and 
already-known sources are labelled, and their relationships are 
discussed in the text. 



